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Abstract 

Spitzer  Infrared  Spectrograph  observations  of  the  Deep  Impact  experiment  in  July  2005  have  created  a  new  paradigm  for  understanding  the 
infrared  spectroscopy  of  primitive  solar  nebular  (PSN)  material — the  ejecta  spectrum  is  the  most  detailed  ever  observed  in  cometary  material. 
Here  we  take  the  composition  model  for  the  material  excavated  from  Comet  9P/Tempel  l’s  interior  and  successfully  apply  it  to  Infrared  Space 
Observatory  spectra  of  material  emitted  from  Comet  C/1995  Ol  (Hale-Bopp)  and  the  circumstellar  material  found  around  the  young  stellar 
object  HD  100546.  Comparison  of  our  results  with  analyses  of  the  cometary  material  returned  by  the  Stardust  spacecraft  from  Comet  81P/Wild 
2,  the  in  situ  Halley  flyby  measurements,  and  the  Deep  Impact  data  return  provides  a  fundamental  cross-check  for  the  spectral  decomposition 
model 8  presented  here.  We  find  similar  emission  signatures  due  to  silicates,  carbonates,  phyllosilicates,  water  ice,  amorphous  carbon,  and  sulfides 
in  the  two  ISO-observed  systems  but  there  are  significant  differences  as  well.  Compared  to  Tempel  1,  no  Fe-rich  olivines  and  few  crystalline 
pyroxenes  are  found  in  Hale-Bopp  and  HD  100546.  The  YSO  also  lacks  amorphous  olivine,  while  being  super- rich  in  amorphous  pyroxene.  All 
three  systems  show  substantial  emission  due  to  polycyclic  aromatic  hydrocarbons.  The  silicate  and  PAH  material  in  Hale-Bopp  is  clearly  less 
processed  than  in  Tempel  1,  indicating  an  earlier  age  of  formation  for  Hale-Bopp.  The  observed  material  around  HD  1 00546  is  located  ^13  AU 
from  the  central  source,  and  demonstrates  an  unusual  composition  due  to  either  a  very  different,  non-solar  starting  mix  of  silicates  or  due  to  disk 
material  processing  during  formation  of  the  interior  disk  cavity  and  planet(s)  in  the  system. 

©  2006  Elsevier  Inc.  All  rights  reserved. 

Keywords:  Infrared  observations;  Spectroscopy  ;  Comets,  origin;  Comets,  composition;  Solar  nebula 


1.  Introduction 

The  recent  Deep  Impact  (DI)  mission  to  the  Comet  9P/ 
Tempel  1  (Tl)  allowed  unprecedented  observations  of  the  pris 
tine  nuclear  material  hidden  beneath  the  processed  surface  layer 
in  typical  comets.  On  2005  July  4,  a  364  kg  impactor  spacecraft 
from  DI  hit  Tl  at  a  relative  velocity  of  10.2  km  s-1 ,  excavating 
~106  kg  of  material  that  was  observed  by  the  instruments  on 
the  DI  flyby  spacecraft  (A’ Hearn  et  al.,  2005),  as  well  as  in  the 
mid-infrared  (TR)  by  the  Spitzer  Space  Telescope  (Werner  et  al.. 
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2004;  Lisse  et  al.,  2006).  The  presence  of  abundant  water  ice  in 
the  ejecta  (Sunshine  et  al.,  2006),  the  minimal  temperature  ex¬ 
cursion  seen  for  the  excavated  dust  (Groussin  et  al.,  2007),  and 
the  arguments  of  Melosh  (1989)  on  impact  cratering  indicate 
that  the  bulk  of  material  released  from  Tl  due  to  DI  was  gently 
removed.  It  may  in  fact  have  been  emitted  in  a  process  simi¬ 
lar  in  intensity  to  emission  mechanisms  observed  in  Comet  C/ 
1995  Ol  (Hale-Bopp).  Fig.  1  compares  the  mid-IR  spectrum  of 
the  ejecta  to  that  of  the  ambient  coma  (Lisse  et  al.,  2006)  taken 
with  Spitzer’ s  Infrared  Spectrograph  (IRS;  Houck  et  al.,  2004) 
and  to  that  of  Hale-Bopp  (HB)  taken  with  the  Short  Wavelength 
Spectrometer  (SWS;  Leech  et  al.,  2003)  on  the  Infrared  Space 
Observatory  (TSO;  Kessler  et  al.,  2003).  Comparison  of  the  two 


70 


C.M.  Lisse  etaL  /  Icarus  187  (2007)  60-86 


Fig.  1 .  Sparer  5-35  pm  spectra  of  Tcmpel  1.  From  bottom  to  top,  spectrum  of  the  ambient  coma,  taken  23  h  before  impact;  spectrum  of  the  ejecta  at  i  —  0.67  h 
after  impact;  ISO  spectrum  of  CVHale-Bopp  (after  Crovisier  et  a  1„  1997);  ISO  spectrum  of  YSO  HD 100546  (after  Malfait  et  al..  1998).  The  spectra  have  been 
divided  by  factor  of  10*  1. 2.7*  and  10  for  ease  of  presentation.  Note  the  logarithmic  scale,  and  the  sharper  resolution  of  features  in  the  Tempel  I  ejecta  spectrum 
vs  the  Hale-Bopp  and  HD  100546  spectra,  h  is  this  increased  feature  resolution  -  Spitzer’s  high  spectral  sensitivity  (compare  the  error  bars  in  the  anissivity 
plots  of  Figs.  2-4)  that  allows  the  much  improved  compositional  determination  of  the  material  in  Comet  9P/Tempel  1  as  compared  to  any  other  comet,  including 
C/Hale-Bopp  (after  Lisse  et  al.,  2006).  The  discontinuities  in  the  IRS  spectra  at  7 5  and  20.5  fan  are  iastiumental  artifacts  due  to  slight  mismatches  in  the  spectra 
from  the  different  IRS  modules. 


comets  shows  numerous  points  of  similarity,  although  the  T1 
spectrum  has  more  structure  in  its  spectrum  than  HB. 

The  1995-1996  apparition  of  Comet  HB,  with  its  strong  out¬ 
flow  of  small,  solid  dust  (Williams  et  al.,  1997;  Lisse  et  al., 
1999;  Harker  et  al.,  2002),  provided  a  “natural”  IR  spectroscopy 
experiment  for  ISO  to  observe.  The  amount  of  large  particulate 
material  (>10  pm)  was  more  abundant  than  in  the  T1  ejecta, 
leading  to  a  smaller  emission  feature  to  continuum  ratio  in  the 
HB  spectra.  On  the  other  hand,  there  was  an  abnormally  large 
amount  of  micrometer-  and  sub  micrometer- sized  solid  mater¬ 
ial  in  the  coma  to  provide  a  highly  featured  spectrum  (Fig.  1 ; 
Crovisier  et  al.,  1997;  Lellouch  et  al.,  1998).  The  presence  of 
this  incredibly  small  and  solid  dust,  in  light  of  the  very  high  rate 
of  emission  of  material  [^lO30  mol  s  1  water  at  1  AU  heliocen¬ 
tric  distance  (Biver  et  al.,  1999, 2002;  Dello  Russo  et  al.,  2000); 
104  kgs-1  dust  (Li  and  Greenberg,  1998;  Lisse  et  al.,  1999; 
Jewitt  and  Matthews,  1999)],  begs  the  question  of  whether  the 
type  of  dust  normally  emitted  by  comets  (Huebner  and  Benkoff, 
1999;  Lisse  et  al.,  2002)  was  altered  by  the  violence  of  the  out¬ 
flow.  Alternatively,  the  dust  from  HB  may  have  been  unusually 
volatile  rich  and  fragile,  although  dust-togas  ratio  estimates  for 
the  comet  are  some  of  the  highest  recorded,  ~3-7.  Toth  and 
Lisse  (2006)  argue  that  Hale-Bopp  is  rotational ly  unstable,  with 
centrifugal  forces  exceeding  the  shear  strength  of  the  nuclear 
material,  leading  to  large  deep  cracks  in  the  body  and  exposure 
of  deep,  volatile-rich  interior  material  to  solar  insolation. 

None  of  the  previously  published  models  for  the  solid  ma¬ 
terial  emitted  by  HB  fit  the  IRS  observations  of  the  material 


excavated  from  T1  at  the  95%  confidence  level  (C.L.;  Lisse 
et  al.,  2006).  In  this  paper,  we  briefly  recap  the  new  compo¬ 
sitional  model  from  the  Dl-Spitzer  observations,  focusing  on 
the  sanity  checks  and  correlations  applied  to  test  the  physical 
reasonableness  of  the  results.  We  then  revisit  the  ISO  spectrum 
of  HB  using  the  new  model.  We  find  that  the  spectra  of  both 
comets  can  be  well  fit  using  a  compositional  model  with  mini¬ 
mal  assumptions  based  on  thermal  emissivity  measurements  of 
powdered  mineral  specimens. 

The  next  step  is  to  extend  the  model  to  other  systems  dom 
mated  by  comet-like  material.  Here,  we  present  a  model  for 
the  ISO  spectrum  of  HD  100546  (HD),  a  ~10  Myr  old  Her 
big  Be9  V  star  which  has  often  been  compared  to  HB  and  other 
comets  (e.g.,  Waelkens  et  al  .,  1 997).  The  age  of  HD  corresponds 
to  the  era  in  which  planetesimals  and  planets  were  formed  in  the 
Solar  System,  and  recent  work  by  Liu  et  al.  (2003)  and  Grady 
et  al.  (2005)  has  demonstrated  disk  clearing  inside  13  AU,  most 
likely  due  to  planet  formation.  We  can  easily  fit  the  observed 
HD  spectrum  with  the  Dl-Spitzer  model  as  well. 

In  this  paper,  we  show  that  the  three  objects  have  similar 
IR  spectra  due  to  the  presence  of  common  fundamental  build 
ing  blocks  of  planetary  systems  (silicates,  water,  amorphous 
carbon,  carbonates,  sulfides,  PAHs),  but  that  under  closer  ex¬ 
amination  they  also  have  significant  differences.  We  then  go  on 
to  show  that  HD  resembles  HB  more  than  T1  in  its  dust  com 
position,  mean  temperature,  and  size  distribution.  From  this, 
we  deduce  that  the  time  of  formation,  more  than  location  of 
formation,  has  the  strongest  effect  on  the  observed  differences 
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and  put  T1  into  context  as  a  body  formed  relatively  late  in  the 
history  of  the  early  Solar  System.  Finally,  we  speculate  from 
the  observed  compositional  trends  on  possible  new  avenues  of 
research  to  be  followed  in  comparative  mineralogy  when  mod* 
eling  the  plethora  of  new  Spitzer  IRS  comet  and  exo- system 
spectra. 

2.  Deep  Impact  and  the  Tfempel  1  dust  model 

The  Spitzer  IRS  observed  T1  on  2-24  July  2005.  The  5- 
35  pm  spectra  were  taken  prior  to,  during,  and  after  the  Deep 
Impact  encounter,  which  occurred  when  T1  was  1.51  AU  from 
the  Sun.  The  material  that  was  ejected  from  depths  as  large  as 
30  m  was  pristine  and  largely  unaltered,  due  to  the  weakness  of 
the  cometary  material  and  the  low  escape  velocity  (~1  ms"1 2 3 4) 
from  the  cometary  surface  (A’ Hearn  et  al.,  2005).  At  the  same 
time,  it  was  de-aggregated  from  loosely  held  fractal  particles 
into  individual  sub-fractal  components  (A’ Hearn  et  al.,  2005; 
Lisse  et  al.,  2006;  Sunshine  et  al.,  2006).  The  observed  material 
had  cooled  from  effects  due  to  the  impact  within  minutes,  and 
the  separation  of  the  ejecta  from  the  ambient  coma  dust  was 
cleanly  made.  The  resulting  highly  structured  spectrum  of  the 
ejecta  shows  over  16  distinct  spectral  features  at  flux  levels  of  a 
few  Jansky s  (Fig.  1)  that  persisted  for  more  than  20  h  after  the 
impact. 

Tb  determine  the  mineral  composition  of  T1  from  the  IRS 
spectra,  we  applied  a  model  based  on  measured  thermal  in¬ 
frared  emission  spectra  for  micrometer- sized  particles  of  over 
80  mineral  species.  Details  of  the  dust  model  are  described  in 
the  Supplementary  online  material  for  Lisse  et  al.  (2006).  This 
model  has  limited  input  assumptions  while  simultaneously  min¬ 
imizing  the  number  of  model  parameters.  Model  solutions  were 
found  by  a  direct  search  through  (composition,  temperature, 
size  distribution,  porosity)  phase  space.  Acceptable  models  for 
the  IRS  spectra  were  required  to  have  Xv  <  1-13  for  325  de¬ 
grees  of  freedom  at  the  95%  Cl.  The  phase  space  search  ex¬ 
cluded  the  vast  majority  of  mineral  species  from  the  T1  ejecta. 
Emission  signatures  due  to  Ca/Fe/Mg  rich  silicates,  carbonates, 
phyllosilicates,  water  ice,  amoiphous  carbon,  and  sulfides  were 
definitively  found  (Fig.  2  and  Table  1 ;  Lisse  et  al.,  2006). 

To  ensure  the  robustness  of  the  T1  ejecta  modeling,  a  number 
of  validation  checks  were  applied: 

1.  All  detected  materials  are  associated  with  major  species 
found  in  the  Interplanetary  Dust  (IPD)  observational  data¬ 
bases  (Bradley,  2002,  2003;  Horz  et  al.,  2000)  and  in  pre¬ 
vious  observations  of  cometary  and  YSO  dust  (Bregman 
et  al.,  1987;  Hanner  et  al.,  1994,  1998;  Wooden  et  al., 
1999,  2000). 

2.  The  elemental  abundances  derived  for  the  observed  refrac¬ 
tory  material  are  consistent  with  Solar  System  abundances. 

3.  The  resulting  particle  size  distribution  agrees  closely  with 
measurements  by  the  DI  instruments  (Sugita  et  al.,  2005; 
Harker  et  al.,  2005;  Sunshine  et  al.,  2006). 

4.  The  derived  dust  temperature  is  consistent  with  the  tem¬ 
perature  measured  by  the  DI  High  Resolution  Instrument 
(HRI)  IR  spectrometer  (Groussin  et  al.,  2007)  and  by 


the  NASA  Infrared  Telescope  Facility  SpeX  instrument 
(Fernandez  et  al.,  2007). 

5.  The  T1  model  results  are  consistent  with  the  Halley  flyby 
evidence  for  silicates,  water,  sulfides,  and  carbonates  to  the 
limit  of  their  sensitivity  (cf.  Clark  et  al.,  1 987;  Jessberger  et 
al.,  1988),  as  discussed  in  Lisse  et  al.  (2006).  The  CHON 
particles  detected  in  the  Halley  flybys  (~25%  by  number) 
are  identified  with  the  dark  amorphous  carbon  component 
(~10%  by  surface  area  in  the  T1  ejecta)  used  here. 

6.  Comparison  of  the  compositional  results  to  the  preliminary 
mineralogy  results  from  the  Stardust  mission  (cf.  Brown 
lee  et  al.,  2006a,  2006b;  Joswiak  and  Brownlee,  2006; 
Zolensky  et  al.,  2006)  finds  consistency  with  the  refrac¬ 
tory  species  in  the  T1  model,  at  the  gross  level.  In  the 
small  number  of  Stardust  terminal  particles  examined  to 
date,  evidence  for  Mg-  and  Fe-rich  olivines,  as  well  as 
Mg-,  Fe-,  and  Ca-rich  pyroxenes  were  found.  Glassy  sili¬ 
cates  and  carbon  were  abundant,  as  were  metal  sulfides.  (It 
must  be  noted  that  there  are  some  differences  between  the 
T1  model  and  the  Stardust  initial  findings  for  the  simple 
terminal  particles,  as  there  is  no  evidence  to  date  of  car¬ 
bonates  or  phyllosilicates  in  the  preliminary  examinations 
of  the  samples  returned  from  Wild  2,  and  the  metal  sul¬ 
fides  appear  most  like  a  non-terrestrial  mixture  of  pyrrhotite 
(Fei  _,S)  and  pendandite  (FeNiS).  On  the  other  hand,  the 
Stardust  initial  analyses  have  shown  a  high  degree  of  parti¬ 
cle  to  particle  variation,  and  a  large  fraction  of  the  deposited 
material  is  found  heated,  altered  and  intermixed  with  the 
silicate  based  aerogel  capture  medium,  complicating  the 
analysis  of  the  bulk  of  the  deposited  material.  Careful  study 
of  the  systematic  effects  in  the  sample  return  should  help 
solve  many  of  these  problems  of  interpretation,  but  it  is 
not  clear  at  the  time  of  this  writing  if  issues  of  small  num¬ 
ber  statistics,  hot  sample  collection  into  aerogel,  collection 
of  mainly  surface  mantle  (and  thereby  altered)  material,  or 
samples  from  a  substantially  younger  comet  are  important 
for  the  comparison.  [For  more  discussion  on  these  points, 
see  Section  5.4.]) 

The  new  modeling  presented  here  allows  us  to  get  beyond 
the  classical  olivine-pyroxene-amorphous  carbon  composition 
to  the  second-order,  less  emissive  species  like  water,  sulfides, 
PAHs,  phyllosilicates,  and  carbonates.  On  the  other  hand,  there 
are  limitations  to  the  method.  There  is  no  petrological  or  iso¬ 
topic  information,  and  the  results  returned  are  for  bulk  aver¬ 
ages  of  the  observed  systems.  For  example,  only  very  abundant 
species  with  strong  emission  features  (>10%  of  the  silicate 
emission  peaks)  will  be  detectable.  Alumina  (AI2O3)  and  hi- 
bonite  (CaAhC^),  extremely  stable  oxides  of  aluminum,  were 
not  definitively  detected  in  the  T1  ejecta  even  though  they 
are  known  minor  constituents  of  the  pre  solar  grain  population 
found  in  interplanetary  dust  particles  (IDPs).  Also,  in  order  to 
cover  the  large  phase  space  of  possible  minerals  present,  we 
assume  a  linear  mix  of  the  extreme  endstates  of  the  mineral 
each  mineral  system,  and  a  linear  shift  in  the  band  positions 
and  strengths  between  the  endstates;  e.g.,  we  linearly  adjust  the 
balance  of  forsterite  (FolOO,  or  MgSi2C>4)  and  fayalite  (FolOO, 
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Fig.  2.  Emissivity  spectrum  for  the  Tempel  1  ejecta,  at  impact  +  45  min,  showing  the  individual  components  of  the  dust  that  create  the  observed  emission,  (a)  To¬ 
tal  spectrum.  The  comet  was  0.75  AU  from  Spitzer  and  1.51  AU  from  the  Sun  at  the  time  of  observation.  The  vertical  scale  is  linear  to  emphasize  the  closeness 
of  the  model  fit  Black;  SST  ejecta  spectrum  divided  by  the  pre  impact  spectrum  Error  bars  are  ~2 a.  Orange  dashed  line  best  fit  model  spectrum.  Colored 
curves — emission  spectra  for  the  constituent  species,  scaled  by  the  ratio  #x(^dusi(a)/)/ ^x(^post-impuct).  with  ^poxt- impact  —  390  K,  and  offset  for  clanty.  Pur¬ 
ple— amorphous  silicates  of  pyroxene  or  olivine  composition.  Light  blue— crystalline  pyroxenes:  ferrosilite,  diopside,  and  orthoenstatite,  in  order  of  20  pm 
amplitude.  Dark  blue-crystalline  olivines,  forsterite  and  feyalite,  in  order  of  20  pm  amplitude.  Red— amorphous  carhon.  Deep  orange— water  ice.  Light  or¬ 
ange  water  gas.  Yellow  PAHs.  Bright  green  carbonates:  siderite  and  magnesite,  by  order  of  7  pm  emissivity  amplitude.  Olive  green  sulfides,  represented  here 
by  niningerite.  Blue  phyllosilicates  (near  the  sulfides),  represented  here  by  the  smectite  nontronite.  (b)  Spectrum  after  subtraction  of  the  best-fit  silicate  compo 
nents.  The  non-silicate  species  best-fit  emission  curves  have  been  overlaid,  scaled  by  a  factor  of  2.  The  sharp  residual  at  19  20  pm  is  due  to  a  poor  match  of  our 
forsterite  emission  spectrum  to  the  data.  We  see  this  mismatch  for  all  three  systems.  (Both  figures  copyright  Lisse  et  aL  2006.) 
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Table  1 

Composition  of  the  best  fit  model  to  the  SST  IRS  Tempel  1  ejecta  spectrum 


Species 

Weighted 
surface  area 

Density 

(gcm~3) 

M.W. 

A'molea 

(rel.) 

Model  T,™ 

(K) 

Model  xi 
if  not  included 

Olivines 

Arnorph  olivine  (MgFe$i04) 

0.17 

3.6 

172 

0.35 

340 

5.92 

Forsterite  (Mg2Si04) 

0.31 

3  2 

140 

0.70 

340 

4.28 

Fayalite  (Fe2Si04) 

0.086 

4.3 

204 

0.18 

340 

1.40 

Pyroxenes 

Arnorph  pyroxene  (MgFeS^Ofc) 

0.041 

3.5 

232 

0.06 

340 

1.38 

Ferrosilite  (F^SijO^) 

0.33 

4.0 

264 

0.50 

295 

9.30 

Diopside  (CaMgS^O*,) 

0.115 

3.3 

216 

0.18 

340 

1.86 

Orthoenstatite  (N^S^O^) 

0.10 

3.2 

200 

0.16 

340 

1.70 

Phyllosilicates 

Smectite  nontronite 

0.14 

2.3 

496 

0.07 

340 

3.76 

Nao.33  Fe2(SUAI)4O|0(OH)2  3H2O 

Carbonates 

Magnesite  (MgC03) 

0.030 

3.1 

84 

0.11 

340 

1.30 

Siderite  (FeCC>3) 

0.051 

3.9 

116 

0.17 

340 

1.83 

Metal  sulfides 

Niningerite  (MgioFegoS) 

0.15 

4.5 

84 

0.92 

340 

2.51 

Water 

Water  ice  (H2O) 

0.049 

1.0 

18 

0.27 

220 

1.40 

Water  ga-s  (H2O) 

0.028 

1.0 

18 

23.7 

220 

1.22 

'Organics ' 

Arnorph  carbon  (C) 

0.068 

2.5 

12 

1.45 

390 

13.3 

PAH(C,oH,4) 

0.039 

1.0 

(178) 

0.022 

N/A 

158 

Note..  ^nw>lc»(0  ~  density (/)/mo»ecular  weight# )  x  surface  area  weighting^);  errors  on  the  abundances  (2a).  Total  best-fit  model  xi  =  1-06.  1TR  at 

1.51  AU  =  230  K_  Table  from  Lisse  et  al.  (2006). 


or  Fe2Si2C>4)  to  fit  the  observed  spectrum,  allowing  us  to  deter¬ 
mine  the  total  number  of  each  atom  present,  but  cannot  distin¬ 
guish  between  the  presence  of  Fo50  (FeMgSi2C>4)  and  a  50-50 
mix  of  (FolOO  +  FalOO).  The  values  in  the  compositional  tables 
should  be  interpreted  in  this  way.  We  also  cannot  distinguish 
easily  between  “glassy  silicate  of  non  stoichiometric  but  near 
olivine  (or  pyroxene)  composition”  and  “amorphous  silicate 
of  olivine  (or  pyroxene)  composition”  and  so  we  assume  the 
presence  of  stoichiometric  glasses  when  modeling  the  glassy 
silicates. 

Despite  these  limitations,  at  this  point  in  time  we  are  able 
to  determine  the  overall  amounts  of  the  different  classes  of  ma¬ 
terial  (olivines,  pyroxenes,  sulfides,  water,  etc.)  and  the  bulk 
elemental  abundances  for  the  most  abundant  atoms.  Applying 
the  model,  with  strong  checks  of  its  validity,  has  great  potential 
for  interpreting  new  mid  IR  spectra  of  distant  dusty  systems 
like  YSOs,  debris  disks,  and  PNs.  We  also  expect  to  update  and 
improve  the  compositional  model  presented  here,  with  its  linear 
sum  of  pure  mineral  end  members,  by  testing,  vs  the  T1  ejecta 
emission  spectrum,  a  more  precise  mineralogical  and  chemical 
formulation  based  on  aggregate  Stardust  measurements  once 
they  are  available  in  the  next  3  to  5  years. 

3.  Modeling  results  for  Comet  Hale-Bopp 

The  2.4-45  pm  ISO  spectrum  of  HB  obtained  by  Crovisier 
et  al.  (1997,  2000)  on  7  October  1996  (heliocentric  distance 


r  =  2.8  AU,  geocentric  distance  A  =  3.0  AU)  has  previously 
been  modeled  by  others  (e.g.,  Malfait  et  al.,  1998;  Wooden  et 
al.,  1999,  2000;  Min  et  al.,  2005).  In  those  analyses,  though, 
the  authors  typically  focus  on  a  specific  species  or  lack  defini¬ 
tive  measures  of  the  basic  nature  and  makeup  of  the  observed 
dust  A  mix  of  material  consisting  of  olivines,  amorphous  car¬ 
bon,  native  iron,  and/or  water  ice  is  usually  assumed,  often 
without  regard  to  cosmic  or  Solar  System  abundances.  Two 
separate  particle  size  distributions  are  often  used  for  the  con¬ 
tinuum  and  spectral  features,  and  a  run  of  temperature  and 
density  is  found  as  a  solution.  Inferences  are  then  drawn  con¬ 
cerning  the  dominant  physical  processes  in  the  system  from 
the  observed  distributions.  Spectral  features  which  cannot  be 
fit,  including  emissions  from  carbonates,  PAHs,  sulfides  and 
phyllosilicates,  are  left  as  a  residual,  and  rarely  if  ever  is  the 
chi-squared  statistic  for  the  presented  best-fit  model  given.  Here 
we  revisit  this  modeling,  using  the  new  results  on  cometary 
composition  found  by  modeling  the  T1  ejecta,  and  a  statisti¬ 
cally  robust  fitting  method.  The  higher  feature-to-continuum 
ratio  and  sensitivity  in  the  Spitzer  spectrum  allowed  for  the 
detection  and  identification  of  more  materials  than  were  previ¬ 
ously  found  in  the  HB  data  alone.  This  accurate  determination 
of  the  composition  of  the  T1  material  enables  us  to  revisit  the 
HB  spectrum  and  remove  some  of  the  ambiguities  produced 
by  the  presence  of  large,  optically  thick  dust  grains  in  the  HB 
coma. 
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Fig.  3.  (a)  Emits  ivity  spectrum  for  the  Hale-Bopp  naturally  emitted  coma,  as  measured  at  2.8  AO  from  the  Earth  and  ISO.  Hiror  bars  are  —2a.  Black— ISO  coma 
spectrum,  divided  by  a  210  K  blackbody.  The  colored  curves  have  the  same  sense  as  in  Fig.  2.  Except  for  Fe-rich  olivine  and  water  gas,  the  entire  suite  of  species 
found  in  the  Tempel  1  ejecta  is  also  found  in  the  Hale  Bopp  coma,  (b)  Spectrum  after  subtraction  of  the  best  fit  silicate  components,  as  in  (a). 


3.1.  Size  distribution  and  mass 

The  new  modeling  was  greatly  aided  by  the  determination 
that  the  particle  size  distribution  of  the  HB  coma  is  dominated 
by  small,  solid,  icy  grains  with  a  power  law  size  distribution 


d n/da~a  3  6  lo4*°  (cf.  Williams  et  al.,  1997;  Lisse  et  al.,  1999; 
Harker  et  al.,  2002),  and  porosity  M).  Using  this  PSD  range, 
we  find  the  total  mass  associated  with  the  spectrum  is  approx 
imately  2  x  109  kg,  three  orders  of  magnitude  more  than  was 
observed  in  the  T1  ejecta.  Even  so,  this  is  still  small  compared 
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with  the  amount  of  dust  seen  in  the  whole  coma  4  tail  structure 
of  HB  (-10n  kg;  Lisse  et  al.,  1998;  Lisse,  2002).  The  differ 
ence  is  due  to  two  effects:  the  small  angular  size  of  the  ISO 
SWS  spectrometer  beam  versus  the  very  extensive  HB  coma  + 
tail;  and  the  fact  that  the  HB  coma,  unlike  the  comae  of  almost 
any  other  comet,  was  optically  thick  near  the  nucleus  due  to  its 
huge  rate  of  emission  of  material. 

3.2.  Temperature 

The  best-fit  simple  blackbody  to  the  continuum  has  a  tem¬ 
perature  of  210  K.  From  our  more  detailed  T1  compositional 
model,  we  find  the  temperature  of  the  smallest  dust  particles 
(0.1- 1.0  pm),  which  superheat  above  local  thermal  equilibrium 
(LTE),  is  —220  K  for  the  olivines,  1 90  K  for  the  pyroxenes,  and 
410  K  for  the  amorphous  carbon.  The  largest  particles  included 
in  the  model,  1000  pm,  are  set  to  7lte  =  282/V(2.8  AU)  = 
169  K,  and  dust  of  intermediate  size  is  scaled  between  the  two 
extremes  (Lisse  et  al.,  2006).  We  find  the  best-fit  temperature  of 
the  water  ice  and  gas,  cooled  by  the  heat  of  sublimation  below 
LTE,  to  be -150  K,  regardless  of  particle  size. 

We  note  here  an  important  point  in  analyzing  the  data:  the 
scattered  sunlight  from  the  coma  of  HB  must  be  carefully  fit 
and  removed  before  modeling  the  data.  Otherwise  the  comet 
appears  to  have  a  huge  super  thermal  continuum.  For  the  large 
majority  of  comets,  the  scattered  light  component  would  be 
negligible  longwards  of  —3  pm.  However,  the  albedo  of  HB’s 
dust  was  so  high  (—40%  vs  —8%  for  Tl’s  ambient  coma)  and 
the  comet  was  observed  relatively  far  from  the  Sun  (—2.8  AU 
vs  a  more  typical  0.7-1. 5  AU)  with  cold  coma  dust,  that  the 
scattered  light  component  was  significantly  enhanced  and  an 
important  component  of  the  measured  5-10  pm  radiation. 

3.3.  Composition 

With  compositional  phase  space  for  comet- like  material  now 
defined  robustly  by  the  DI  results,  and  porosity  and  size  distrib¬ 
ution  given  by  the  literature,  we  solve  iteratively  for  the  relative 
amount  of  each  species  while  demanding  good  agreement  with 
the  measured  ISO  spectrum.  We  find  that  there  is  no  need  to 
make  the  ad  hoc  assumption  of  two  separate  particle  size  dis¬ 
tributions,  hollow  spheres,  long  needles,  or  an  unusual  mineral 
composition.  The  same  material  that  yields  the  sharp  spectral 
features  also  provides  the  continuum,  from  its  larger,  optically 
thick  pieces.  In  fact,  the  number  of  materials  making  up  the  ma¬ 
jority  of  the  coma  dust  of  HB  is  fewer  than  for  T1 . 

The  results  of  fitting  the  SWS  spectrum  of  HB  taken  at  2.8 
AU  are  given  in  Table  2,  Fig.  3a  (full  spectrum),  and  Fig.  3b 
(residuals  after  the  dominant  silicate  and  phyllosilicate  best-fit 
model  is  removed).  The  best-fit  model  had  a  x?  =  0.99;  the 
last  column  of  the  table  indicates  the  change  in  xl  if  a  par¬ 
ticular  species  is  not  included  in  the  model.  Fig.  3  shows  the 
ISO  data  with  best-fit  spectral  model  superimposed,  as  well 
as  the  spectrum  of  each  contributing  species  broken  out  sep¬ 
arately.  In  contrast  to  T1  (Figs.  2,  5),  there  is  no  Fe-rich  olivine 
(fayalite,  Fe2Si04).  There  is  also  more  amorphous  olivine  rel¬ 
ative  to  Mg-rich  crystalline  forsterite  (Mg2SK>4),  as  was  found 


by  Wooden  et  al.  (1999).  Total  pyroxenes  (e.g.,  MgFeSi206), 
both  crystalline  and  amorphous,  are  present  in  about  the  same 
amount  as  for  Tl,  consistent  with  the  results  of  Wooden  et  al. 
(1999,  2000).  The  fraction  of  amorphous  pyroxene  material  is 
much  higher  in  HB  than  Tl. 

Phyllosilicates  (as  represented  by  the  smectite  nontron- 
ite,  Nao.33Fe2(Si,Al)40io(OH)2'3H20)  seem  slightly  enhanced 
and  carbonates  (magnesite  MgCC>3  and  siderite  FeCC>3)  ex¬ 
tremely  enhanced,  by  a  factor  of  3,  relative  to  Tl.  The  latter 
statement  depends  strongly  on  the  veracity  of  the  ISO  SWS 
5-8  pm  data  where  the  error  bars  are  large.  However,  the  chi- 
squared  fitting  routine  used  to  derive  the  composition  takes 
these  effects  into  account  and  omitting  the  carbonate  from  the 
model  more  than  doubles  x2  (Table  2).  The  sulfide  emission 
feature  (niningerite  [Mg,  Fe2 1  ]S)  is  similar  in  the  two  comets, 
although  the  best-fit  ningerite  composition  is  Feo.9Mgo.1S  for 
HB  (and  HD,  see  Section  4.2)  rather  than  the  Feo.5Mgo.5S  found 
for  Tl.  The  water  ice  signal  is  much  larger  in  HB  than  for  Tl 
because  HB  was  beyond  the  ice  line  at  the  time  of  observation. 
Water  ice  is  thus  more  stable,  and  the  effective  temperature  we 
find  for  the  water  ice,  —140  K,  is  near  the  LTE  temperature. 

The  detection  of  PAHs  in  the  comet  deserves  some  dis¬ 
cussion.  We  find  that  in  the  best-fit  HB  model,  there  is  about 
twice  as  much  PAH  material  per  mole  of  dust  than  was  found 
for  Tl.  For  such  an  active  comet  with  strong  signatures  of  or¬ 
ganic  species  such  as  methane,  ethane,  and  methanol  (Biver 
et  al.,  1999;  Bocketee-Morvan  et  al.,  2000)  and  a  significant 
amount  of  amorphous  carbon,  the  presence  of  PAHs  was  not 
unexpected.  The  3.3  pm  PAH  feature  was  tentatively  detected 
amongst  OH  and  H2O  lines  in  a  few  comets,  P/Halley,  Levy, 
and  Brorsen-Metcalf  (Encrenaz  et  al.,  1988;  Bockelee-Morvan 
et  al.,  1995),  and  fluorescence  features  attributed  to  aromat¬ 
ics  were  also  tentatively  detected  in  the  near-UV  in  P/Halley 
(Moreels  et  al.,  1994;  Clairemidi  et  al.,  2004).  However,  the 
detection  of  PAHs  in  HB  disagrees  with  previous  determina¬ 
tions  for  this  comet  (Cmvisier  et  al.,  2000).  The  difference  in 
the  current  work  is  the  proper  assignment  of  the  compositional 
materials  making  up  the  coma  dust  and  the  careful  removal  of 
the  dominant  silicate  and  amorphous  carbon  emission  from  the 
spectra,  allowing  detection  of  the  secondary  emitters.  Also,  the 
Tl  results  showed  that  ionized  PAHs  dominated  in  the  ejecta 
(Fig.  2b),  implying  an  ionization  timescale  of  less  than  1000  s, 
so  that  we  can  expect  the  HB  coma  PAHs  to  be  ionized  as 
well.  Detection  of  PAHs  using  the  feature  at  3.29  pm  would 
therefore  also  be  problematic,  as  the  amplitude  of  this  emis¬ 
sion  is  decreased  by  about  a  factor  of  20  compared  to  ionized 
PAHs.  The  3.29  pm  flux  we  estimate  from  our  Tl  model,  0.02 
Jy,  corresponds  to  one- fifth  the  amplitude  or  about  ler  of  the 
most  pronounced  PAH  emission  feature  from  Comet  Levy  in 
Bockelee-Morvan  et  al.  (1 995).  The  PAH  signal  for  HB  is  about 
three  times  as  large,  or  3<x. 

3.4.  Elemental  abundances 

Far  and  away  the  most  abundant  atomic  species  in  the  Prim¬ 
itive  Solar  Nebula  (PSN),  as  in  the  ISM,  was  H,  followed  by 
He.  Heavier  and  more  complex  nuclei  are  much  rarer,  by  orders 
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Tabic  2 

Composition  of  the  best-fit  model  to  the  ISO  SWS  2.8  AU  Hale  Bopp  spectmm 


Species 

Weighted 
surface  area 

Density 

(gem-3) 

M.W. 

Wmoles 

Get) 

Model  Traax 
(K) 

Model  xi 

if  not  included 

Olivines 

Amorph  olivine  (MgFeSi04) 

0.18 

3.6 

172 

0.37 

220 

8.76 

Forsterite  (Mg2Si04) 

0.32 

3.2 

140 

0.73 

220 

8.17 

Fayalite  (Fe2S  i04) 

0.03 

4.3 

204 

0.06 

220 

1.07 

Pyroxenes 

Amorph  pyroxene  (MgFeS^Oft) 

0.13 

3.5 

232 

0.20 

220 

4.72 

Ferrosilite  (F^S^Oft) 

0.11 

4.0 

264 

0.17 

190 

3.39 

Diopside  (CaMgSi20^) 

0.05 

3.3 

216 

0.08 

220 

1.22 

Orthoenstatite  (Mg2Si2Of>) 

0.12 

3.2 

200 

0.19 

220 

3.43 

PkyllosUicates 

Smectite  nontranite 

0.18 

2.3 

496 

0.08 

220 

8.15 

Naa33Fe2(SiwM)4O|0(OHh  3H20 

Carbonates 

Magnesite  (MgCOy) 

0.15 

3.1 

84 

0.55 

220 

1.92 

Siderite  (FeCO^) 

0.18 

3.9 

116 

0.60 

220 

2.20 

Metal  sulfides 

Niningerite  (MgjoFe^S) 

0.15 

4.5 

84 

0.80 

220 

4.95 

Water 

Water  ice  (H20) 

0.34 

1.0 

18 

1.89 

150 

8.68 

Water  gas  (H20) 

0.00 

1.0 

18 

0.00 

150 

0.99 

‘Organics’ 

Amorph  caibon  (C) 

0.08 

2.5 

12 

1.67 

410 

1.03 

PAH(C,0Hi4) 

0.14 

1.0 

(178) 

0.08 

N/A 

1.50 

Nat*.  NtwAtssd)  ~  density  (i )/ molecular  weight (i)  x  surface  area  weightmg(i);  =20%  errors  nn  the  abundances  (2cr).  Total  best-fit  model  xl  =0.99.  fXR  at 


2.8  AU  *  169  K. 


of  magnitude,  with  the  abundance  failing  rapidly  with  atomic 
number  Z.  Simply  put,  this  is  a  statistical  effect — the  more  pro¬ 
tons  and  neutrons  required  to  make  up  a  nucleus,  the  lower 
the  probability  the  nucleus  was  formed.  There  are  exceptions 
to  this  rule — the  metals  CNO  are  enhanced  due  to  their  for¬ 
mation  in  the  catalytic  CNO  cycle,  and  Fte  is  super-abundant 
because  it  lies  at  the  minimum  of  the  nuclear  energy  binding 
curve.  Similarly,  nuclei  with  especially  stable  binding  energy 
configurations— like  the  even-even  nuclei  with  half-filled  or 
filled  nuclear  shells — are  extra-abundant. 

The  most  abundant  atomic  species  found  in  primitive  solid 
state  materials  are  H:C:0:Si:Mg:Fe.  Nitrogen,  while  nearly  as 
abundant  as  C  and  O,  resides  mostly  in  the  icy  volatile  species 
N2  and  to  a  lesser  extent  in  NH3  (e.g.,  Womack  et  al.,  1992). 
The  bulk  of  the  H,  C,  and  O  will  be  in  volatile  organic  ices 
and  in  water  ice,  but  a  fraction  will  be  contained  in  the  refrac¬ 
tory  dust.  Since  these  atomic  species  are  much  more  abundant 
than  the  rock-forming  elements  Si,  Mg,  and  Fe,  they  appear 
in  similar  numbers.  The  next  tier  down  of  abundant  species, 
Ca:Al:Na,  are  also  important  rock- forming  elements,  and  Ca 
and  Al  form  extremely  stable  oxides.  These  are  often  harder  to 
detect  though,  because  of  their  lower  solar  abundance. 

For  the  T1  model,  the  relative  elemental  abundance  in  the 
ejected  dust  are  H:C:O.Si:Mg:Fe:S:Ca:Al  ==  0.42:0.58:3.9:1.0: 
0.88:0.79:0.29:0.054:^0.085  (for  Si=  1.0).  For  reference,  the 
solar  abundance  values  vs  Si  are  3.1  x  14.1: 1.00: 1.05: 

0.87:0.43:0.063:0.072  (Anders  and  Grevesse,  1989;  Asplundet 


al.,  2005).  The  Si:Mg:Fe:Ca:Al  ratios  we  find  for  T1  are  con¬ 
sistent  with  Solar  System  abundances  determined  from  the  Sun 
and  from  Cl  chondrites  (Asplund  et  al.,  2005;  Fig.  6).  Hydro 
gen  is  very  deficient  as  expected  for  small  primitive  bodies* 
inability  to  retain  solar  nebula  H2.  Carbon  is  deficient  by  a 
factor  of  ~13.  Oxygen  is  less  deficient  in  the  dust  than  car¬ 
bon,  depleted  by  a  factor  of  ~3.5,  presumably  due  to  its  more 
lithophilic  chemistry.  Sulfur  appears  deficient  by  a  factor  of  1.5, 
which  suggests  that  we  may  not  have  well  characterized  the 
sulfide  constituents  or  that  primitive  bodies  do  not  efficiently 
accrete  the  more  volatile  sulfur  compounds  such  as  H2S  or  SO2 
from  the  nebula. 

For  the  HB  model,  the  relative  elemental  abundance  in  the 
coma  are  H:C:0:Si:Mg:Fe:S:Ca:Al  =  1.6:12:5.3:1.0:1.2:0.86: 
0.34:0.054:  <0.11  (for  Si  =  1.0).  These  abundances  are  gener 
ally  closer  to  the  solar  abundances  (Fig.  6)  than  are  those  forTI . 
More  H  is  present  in  the  form  of  water  ice  for  HB,  as  some  of 
that  volatile  fraction  was  still  in  the  solid  phase  when  the  comet 
was  observed  at  r  =  2.8  AU  where  7lte  =  167  K.  Carbon  is  not 
deficient  at  all.  Oxygen  is  highly  deficient,  with  the  balance  of 
atoms  contained  in  the  large  amount  of  coma  H2O  and  CO/CO2 
gas  in  this  most  active  of  modem  comets.  The  Si:Mg:Fe:Ca:Al 
rock  forming  element  ratios  trends  are  nearly  chondritic.  Sulfur 
appears  deficient  by  approximately  30%,  suggesting  again  that 
we  are  not  well  characterizing  the  sulfide  constituents  or  that 
comets  do  not  efficiently  accrete  volatile  forms  of  S.  A  simi¬ 
lar  suggestion  was  made  by  Lisse  et  al.  (2006)  and  has  been 
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reinforced  by  recent  Stardust  results.  Keller  et  al.  (2006)  and 
Zolensky  et  al.  (2006)  found  pentlandite  rather  than  ningerite 
to  be  the  majority  sulfide  phase  in  the  dust  of  Comet  Wild  2. 

4.  Modeling  results  for  HD  100546 

The  combined  presence  of  crystalline  silicates,  phyllosili- 
cates,  carbonates,  and  icy  volatiles,  and  the  small  amount  of 
PAHs  in  Tl’s  interior  has  fundamental  implications  for  the  na¬ 
ture  of  the  material  incorporated  into  the  comet  4.5  Gyr  ago 
from  the  PSN.  Crystalline  silicates  formed  and  PAHs  were  de¬ 
stroyed  within  the  present  orbit  of  Mercury.  Phyllosilicates  and 
carbonates  formed  in  regions  of  copious  water  between  ~0.5 
and  ~2.5  AU.  Icy  volatiles  were  accreted  as  the  comet  itself 
coalesced  from  these  components  out  beyond  the  giant  planets. 
One  of  the  most  likely  possibilities  is  that  there  was  pronounced 
mixing  in  the  PSN  at  the  time  of  the  comet’s  creation  (Nuth 
et  al.,  2000;  Bocketee-Morvan  et  al.,  2002;  Lisse  et  al.,  2006). 
Alternative  explanations  of  the  observed  chemistry  include  neb¬ 
ular  shockwave  heating  (Harker  and  Desch,  2002)  and  parent 
body  alteration  chemistry  over  the  past  4.5  Gyr  (Lisse  et  al., 
2006). 

With  this  in  mind,  it  is  natural  to  turn  to  mid-IR  observa¬ 
tions  of  other  stellar  systems  thought  to  be  evolving  through 
the  protoplanetary  stage,  dominated  by  cometary  bodies  and 
comet-like  material.  In  the  current  paradigm  of  planetary  sys¬ 
tem  formation,  the  evaporation  and  condensation  of  refractory 
and  icy  materials  takes  place  in  a  thick  protostellar  disk,  fol¬ 
lowed  by  the  settling  of  solid  grains  into  a  thinner  disk  on  a 
timescale  of  105  yr.  Over  the  next  107  yr,  low  velocity  colli¬ 
sions  among  the  grains  result  in  the  accumulation  and  inward 
migration  of  cometary  objects,  with  sizes  ranging  from  mm 
to  km  (Weidenschilling,  1997).  Planetesimal  reservoirs,  i.e., 
protoplanetary  disks,  consistent  with  this  paradigm  have  been 
found  around  a  growing  number  of  young,  T  Tauri  stars  (cf. 
review  by  Dutrey  and  Guilloteau,  2004).  ISO  significantly  in¬ 
creased  our  knowledge  of  these  disks  and  now  Spitzer  is  adding 
new  insight  into  this  phenomenon  (e.g.,  Forrest  et  al.,  2004; 
Megeath  et  al.,  2005). 

The  prototypical  YSO  system  studied  for  its  cometaiy  con¬ 
nections  is  HD  100546,  a  Herbig  Be9  V  star  at  ~10  Myr  age 
since  formation  and  103  pc  distance  (van  den  Ancker  et  al., 
1997).  Herbig  Ae/Be  stars  are  systems  of  intermediate  mass 
still  in  their  pre-main  sequence  (PMS)  phase  of  evolution.  They 
have  lost  most  of  their  surrounding  envelope  of  infalling  gas 
and  dust  but  are  not  yet  fusing  hydrogen  into  helium.  Instead 
they  are  heated  mainly  by  gravitational  contraction  and  material 
infall.  The  ISO  SWS  spectrum  of  this  system  is  often  com¬ 
pared  to  that  of  HB  in  the  literature  (e.g.,  Waelkens  et  al.,  1997; 
Malfait  et  al.,  1998),  and  it  is  used  as  a  laboratory  for  nebular 
formation  processes  because  of  its  apparent  primitive  makeup 
and  young  age  (Grady  et  al.,  2001). 

4. 1.  Temperature,  mass  and  size  distribution 

Here  we  model  the  spectral  decomposition  for  the  ISO  SWS 
spectrum,  as  we  did  for  the  HB  spectrum.  The  major  dif¬ 


ferences  between  the  YSO  and  cometary  cases  are  that  for 
HD  100546  (1)  a  stellar  photospheric  component  must  be  re¬ 
moved,  and  (2)  a  run  of  material  density  and  temperature  vs 
stellar  distance  r*  must  be  addressed.  The  first  issue  is  rela¬ 
tively  straightforward  to  deal  with,  and  we  have  used  a  Ku- 
rucz  model  scaled  to  short  wavelength  photometric  and  2-5  pm 
spectroscopic  measurements  of  the  system  (Hu  et  al.,  1989; 
Waelkens  et  al.,  1997)  to  remove  the  photospheric  contribu¬ 
tion.  If  this  is  not  done  properly,  then  like  the  effect  of  the 
scattered  sunlight  from  HB’s  coma,  the  YSO  appears  to  have 
a  huge  super-thermal  continuum  component. 

To  deal  with  the  second  issue,  we  simplify  the  problem 
by  assuming  that  the  observed  IR  emission  will  be  domi¬ 
nated  by  the  hottest  aod  densest  material.  The  best  simple 
blackbody  fit  to  the  underlying  dust  continuum  appears  to 
come  from  a  hot  region,  at  T  =  250  K,  similar  to  the  ef¬ 
fective  temperature  of  the  zodiacal  dust  as  seen  from  Earth, 
plus  a  cooler  outer  ring  with  effective  temperature  T  =  135  K, 
as  used  in  the  modeling  of  Bouwman  et  al.  (2003).  In  con¬ 
trast,  from  our  detailed  T1  compositional  model,  we  find 
a  large  range  of  dust  particle  temperatures  at  one  astro¬ 
centric  distance,  depending  on  the  particle  composition  and 
size  influencing  its  ability  to  absorb  and  re-radiate  starlight. 
The  best-fit  temperature  of  the  smallest  dust  particles  (0.1- 
1.0  pm)  is  ~250  K  for  the  olivines,  190  K  for  the  pyrox¬ 
enes,  150  K  for  the  water  ice  and  gas,  and  410  K  for  the 
amorphous  carbon.  (Water  ice  and  gas  tend  to  “run  cold” 
compared  to  other  dust  due  to  the  significant  energy  ab¬ 
sorbed  in  vaporizing  water.)  These  results  indicate  an  LTE 
of  about  170  K,  with  the  dominant  emitters,  including  both 
hot  silicaceous  dust  and  cold  icy  dust,  all  co-located  at  the 
same  distance  from  the  primary,  ~ 12-14  AU  (assuming  L*  = 
22Lsoior;  van  den  Ancker  et  al.,  1997;  Weinberger  et  al., 
2003).  We  thus  favor  a  model  of  the  system  with  the  domi¬ 
nant  actively  emitting  material  localized  in  space.  Note  that 
13  AU  is  the  same  distance  as  estimated  by  Grady  et  al. 
(2005)  for  the  inner  wall  of  a  cleared  cavity  in  the  disk  from 
HST/STIS  imaging  observations,  which  would  represent  the 
densest  and  hottest  emitting  dusty  material  in  the  HD  sys¬ 
tem. 

Using  these  dust  temperatures,  the  T1 -based  dust  model 
works  well,  with  Xv  =  1  -07.  The  model  results  are  listed  in  Ta¬ 
ble  3  and  shown  graphically  in  Fig.  4.  The  best- fit  particle  size 
distribution  from  our  modeling,  dn/da  =  a  4  0lt0  05>  fa  steeply 
weighted  towards  the  small,  unaggregated  particle  sizes.  This 
distribution  is  steeper  even  than  the  a~3*8  distribution  found 
for  the  de-aggregated  ejecta  produced  in  the  DI  experiment 
(Lisse  et  al.,  2006)  and  is  steeper  than  all  but  the  most  extreme 
HB  models.  The  total  amount  of  mass  associated  with  the  ISO 
spectrum,  estimated  by  integrating  up  to  the  largest  detectable 
particles  (20  pm  radius),  is  at  least  1022  kg,  about  1  /7  the  mass 
of  the  Earth’s  moon.  Much  larger  total  masses  would  be  esti¬ 
mated  by  assuming  maximum  particle  sizes  on  the  centimeter 
or  meter  scale,  or  if  the  system  is  optically  thick  beyond  the  in¬ 
ner  cavity  wall.  From  both  the  mass  and  temperature  results,  it 
is  clear  that  we  are  observing  a  huge  amount  of  circumstellar 
material  spread  out  over  a  few  AU  in  this  system. 
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Wavelength  (pm) 

(a) 


(b) 

Fig.  4.  Emissivity  spectrum  for  the  ISO  HD  100546  spectrum,  (a)  Total  spectrum.  Color  scheme  is  the  same  as  in  Fig.  2.  Error  bars  are  =2 a.  Black — HD  100546 
excess  spectrum  (i.e.,  central  source  removed),  divided  by  a  250  K  blackbody.  Orange  dashed  line  best-fit  preliminary  model,  fitting  all  major  features  with  the 
1  temperature  assumption.  Unlike  Tempel  1,  the  only  major  olivine  species  found  Is  Mg-ridi  and  crystalline,  while  the  pyroxenes  are  limited  to  Fe-rich  crystalline  and 
the  amorphous  glassy  species.  No  Fe-rich  carbonates  are  detected,  and  the  phyllosilicates  are  relatively  underabundant,  while  PAHs  are  superabundant,  (b)  Spectrum 
after  subtraction  of  best-fit  silicate  components,  as  for  (a). 


4.2.  Composition  spectrum.  The  results  are  given  in  Table  3,  Fig.  4a  (full  spec¬ 

trum),  and  Fig.  4b  (residuals  after  the  dominant  silicate  and 
Applying  the  same  analysis  techniques  as  were  used  for  HB,  phyllosilicate  best-fit  model  is  removed).  What  is  striking  about 

we  obtain  a  good  compositional  fit  to  the  HD  100546  ISO  SWS  the  resulting  overall  HD  100546  composition  is  that  it  appears 
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Table  3 

Composition  of  the  best-fit  model  to  the  ISO  SWS  HD 100546  spectrum 


Species 

Weighted 
surface  area 

Density 

(gem-3) 

M.W. 

^moles 

(rei) 

Model  7niax 
(K) 

Model  xl 
if  not  included 

Olivines 

Araorph  olivine  (MgFeSi04) 

0.00 

3.6 

172 

0.00 

250 

1.07 

Forsterite  (Mg2SiC>4) 

0.35 

3.2 

140 

0.80 

250 

7.36 

Fayalite  (Fe2Si04) 

0.03 

4.3 

204 

0.06 

250 

1.13 

Pyroxenes 

Arooiph  pyroxene  (MgFeS^Ofc) 

0.35 

3.5 

232 

0.53 

250 

20.0 

Ferrosilite  (F^S^O^) 

0.15 

4.0 

264 

0.23 

190 

3.34 

Diopside  (CaMgSi206> 

0.00 

3.3 

216 

0.00 

250 

1.07 

Orthoenstatite  (N^S^Ofc) 

0.03 

3  2 

200 

0.05 

250 

1.16 

Phyllosilicates 

Smectite  nontronite 

0.14 

2.3 

496 

0.07 

250 

4.03 

Nao.33Fe2(SMd)40io(OH)2  3H20 

Carbonates 

Magnesite  (MgCC^) 

0.08 

3.1 

84 

0.30 

250 

1.37 

Siderite  (FeCC^) 

0.00 

3.9 

116 

0.00 

250 

1.07 

Meted  sulfides 

Niningerite  (MgioFe^S) 

0.15 

45 

84 

0.80 

250 

2.31 

Water 

Water  ice  (H20) 

0.46 

1.0 

18 

2.55 

150 

9.66 

Water  gas  (H20) 

0.18 

1.0 

18 

1.007? 

150 

1.14 

4 Organics ' 

Aroorph  carbon  (C) 

0.33 

25 

12 

6.8 

410 

2.51 

PAH  (C]qH]4) 

0.51 

1.0 

(178) 

0.29 

N/A 

4.78 

Note,  tynoic s(i)  ~  density (0/mnlecular  weight(/)  x  surface  area  weightmg(/);  -10%  errors  on  the  abundances  (2o).  Total  best-fit  model  y}  =  107.  I,TR  at 
13  AU  =  170  K. 


in  most  respects  to  be  very  similar  to  that  of  HB,  and  rather  dif¬ 
ferent  than  the  T1  spectrum  (Fig.  5).  Sulfides,  phyllosilicates, 
and  amorphous  carbon  all  seem  to  be  present  in  roughly  the 
same  amounts  as  seen  for  the  comets.  Water  is  present  as  both 
water  ice  and  gas,  indicating  the  observed  material  is  straddling 
the  ice  line  in  the  system.  A  mix  of  water  ice  and  gas  was  seen 
for  Tl,  due  to  the  unusual  circumstances  of  quickly  releasing 
water  ice  at  1 .51  AU  from  the  impact,  but  there  is  abundant  evi¬ 
dence  that  the  ice  began  rapidly  subliming  (e.g.,  Ktippers  et  a!., 
2005;  Sunshine  et  al.,  2006).  For  HB,  only  water  ice  was  found, 
as  expected  for  a  comet  at  2.8  AU  from  the  Sun,  beyond  the 
ice  line.  Compared  to  total  normalized  number  of  moles  of  ma¬ 
terial  found,  the  relative  amount  of  PAHs  in  the  two  systems 
is  found  to  be  the  same  (^2  x  10“2)  and  is  about  a  factor  of 
25  larger  than  the  PAH  abundance  in  Tl.  The  main  difference 
between  HD  and  HB  is  the  super  abundance  of  amorphous  py¬ 
roxene  and  absence  of  amorphous  olivine.  Also,  carbonates  are 
present  as  only  the  Mg-rich  component,  at  about  one-fifth  the 
abundance  seen  in  HB  and  about  one-half  the  abundance  found 
in  Tl. 

The  presence  of  abundant  forsterite,  water  ice,  and  PAHs 
agrees  with  the  findings  of  Malfait  et  al.  (1998)  and  Bouwman 
et  al.  (2003)  but  the  presence  of  sulfides,  crystalline  pyroxenes, 
and  the  lack  of  FeO  do  not.  The  relative  amount  of  phyllosil- 
icate  and  carbonate  is  similar  to  what  is  seen  in  the  comets, 
suggesting  these  species  were  also  present  in  the  PSN. 


4.3.  Elemental  abundances  for  HD  100546 

The  relative  elemental  abundances  derived  for  the  HD 
100546  nebula  are  H:C:0:Si:Mg:Fe:S:Ca:Al  =  3.0:8.4:4. 1:1.0: 
1.1:0.74:0.52:0.01 1:  <0.07  (for  Si  =  1.0).  Compared  to  the  HB 
coma  dust,  which  is  the  closest  of  the  3  systems  to  solar  abun¬ 
dances,  there  is  more  H  present,  in  the  form  of  PAHs.  The  other 
major  H  reservoir,  water  ice,  has  about  the  same  relative  abun¬ 
dance  in  both  systems.  Oxygen  is  markedly  deficient  compared 
to  the  solar  value,  as  for  Tl  and  HB.  In  contrast,  the  C  abun¬ 
dance  is  approximately  that  expected  for  a  solar  abundance  with 
some  carbon  (~1  /3)  found  in  gaseous  species.  Si  and  Mg  are 
close  to  the  solar  ratios;  Fe  is  decreased  in  abundance  by  16%,  a 
value  that  is  within  the  uncertainties  of  the  abundance  estimates. 
The  Ca  abundance  is  markedly  low  vs  solar  and  the  two  comet 
systems.  The  S  component,  as  for  Tl  and  HB,  is  low  by  about 
30%  from  what  would  be  expected  for  a  solar  abundance  star. 

5.  Discussion 

5.1.  Temperature  and  mass  of  the  3  systems 

The  temperature  of  the  Tl  ejecta  was  ~340  K,  as  ex¬ 
pected  for  dust  heated  by  sunlight  at  1 .5  AU,  but  very  unlike 
the  temperatures  (>100K)  in  its  region  of  formation  (r»»  = 
50-200  AU).  The  temperatures  of  the  observed  dust  in  the  disk 
of  HD  100546  (135-250  K)  and  the  coma  of  HB  (~200  K) 
are  very  similar,  close  to  the  sublimation  temperature  of  wa- 
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Fig.  5.  Relative  mineral  abundances  for  the  3  systems.  Tempel  1  is  given  in  orange,  Hale-Bopp  is  denoted  by  blue,  and  HD  100546  by  black.  The  water  (gas  +  ice) 
abundances  have  been  divided  by  a  factor  of  5  and  the  amorphous  carbon  abundances  by  a  factor  of  10  for  ease  of  presentation. 


ter  ice.  The  material  observed  in  HD  100546  corresponds  to 
the  material  that  was  incorporated  into  pianetesimals  in  the 
inner  giant  planet  region  and  into  the  outer  main  belt  aster 
oids.  HB,  an  Oort  cloud  comet,  could  very  well  have  been 
formed  in  a  similar  region  of  the  PSN.  All  3  systems  show 
evidence  for  fine  water  ice  particles,  a  marker  for  “gentle¬ 
ness”  of  the  processing  of  the  observed  material.  None  of 
the  3  systems  was  warm  enough,  nor  did  the  DI  impact  im¬ 
part  sufficient  energy  to  the  cold  ejecta  (A’ Hearn  et  al.,  2005; 
Lisse  et  al.,  2006)  tn  create  any  laTge-scale  changes  in  the  re¬ 
fractory  dust  materials,  although  the  some  of  the  more  volatile 
icy  materials  may  have  been  driven  off,  and  some  of  the  water 
present  vaporized. 

The  observed  dust  masses  of  the  three  objects  vary  widely. 
From  smallest  to  largest,  we  observed  x  106  kg  of  T1  ejecta, 
2  x  109  kg  of  HB  coma,  and  1  x  1022  kg  of  HD  100546  disk  in 
the  0.1-20  um  particle  size  range.  A  number  of  things  become 
clearer  once  this  is  noted.  The  Spitzer  T1  ejecta  measurements 
were  extremely  sensitive,  revealing  a  number  of  new  spectral 
features  despite  observing  103  times  less  material  than  was 
observed  by  ISO  in  HB.  Also,  the  total  amount  of  material  re¬ 
leased  in  the  ejecta  was  only  a  small  fraction,  ~10-4,  of  the 
mass  of  the  dust  in  the  ambient  T1  coma.  Even  so,  the  emission 
from  this  dust  was  a  good  25%  of  the  total  observed  emission 
from  the  comet  while  the  ejecta  remained  near  the  nucleus,  be¬ 
cause  of  the  very  large  emission  efficiency  for  small,  crystalline 
dust  particles.  The  amount  of  HB  material  in  the  ISO  beam  was 
a  relatively  small  fraction  (~1%)  of  the  total  dust  mass  in  the 
coma(~I  x  1011— 1  x  101 2  kg),  but  focused  on  the  freshest  and 
newest  material  emitted  from  the  nucleus.  The  1  x  1022  kg  ob¬ 
served  material  in  HD  100546  is  not  from  a  few  comets,  or  even 
one  giant  comet.  The  amount  of  material  observed  is  equivalent 


to  the  mass  of  ~108  typical  JFC  comets,  or  the  mass  of  a  Kuiper 
belt.  The  system  appears  to  be  a  disk  of  material  with  a  range 
of  particle  sizes  and  the  process  of  planet  formation  going  on 
nearby.  Thus  we  can  expect  the  presence  of  particles  larger  than 
20  um  with  appreciable  mass  but  low  emitting  surface  area  to 
mass  ratio,  so  this  is  most  likely  a  lower  limit  to  the  total  mass 
present. 

5.2.  Composition  of  the  3  systems — overview 

For  an  ~10  Myr  old  system  with  at  least  108  comet’s  worth 
of  mass  spread  out  over  many  AU,  in  the  process  of  disk  clear¬ 
ing  and  planet  building,  it  is  quite  striking  that  HD  resembles 
HB  closely  in  composition  (as  noted  in  the  HD  composition 
section  above,  and  illustrated  in  the  barchart  of  Fig.  5).  As  the 
HD  material  is  located  in  a  region  analogous  to  the  inner  giant 
planet  region  of  the  PSN,  this  would  seem  to  indicate  that  HB 
was  also  formed  in  the  inner  giant  planet  region.  It  is  tempting 
to  “date”  HD  as  younger  than  HB  and  T1  based  on  the  large 
amount  of  amorphous  pyroxenes  and  PAHs,  both  primitive 
species,  and  the  size  distribution  dominated  by  extremely  small 
grains,  implying  little  particle  aggregation.  We  are  given  pause, 
however,  by  the  unusually  small  amount  of  amorphous  olivine, 
the  lack  of  Fe-rich  carbonates,  and  the  evidence  for  nearby 
planet  formation  and  disk  clearing.  HD  is  relatively  well  dated 
as  a  system  of  ~10  Myr  age,  more  evolved  than  other  Herbig 
Ae-Be  stars  by  a  few  Myr  (Grady  et  al.,  2005).  The  giant  planet 
cores  in  the  Solar  System  are  thought  to  have  formed  within  5  to 
1 0  Myr  (Boss,  2004),  so  by  inference,  HB  is  this  age  or  younger. 
The  HD  100546  Be9  V  primary  is  of  very  unusual  composi¬ 
tion,  of  the  “lambda  Bootis”  type  (Acke  and  Waelkens,  2004). 
This  type  of  star  has  the  chemical  peculiarity  of  having  photo - 
spheric  abundances  in  which  C,  N,  O,  and  S  are  roughly  solar 
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Fig.  6.  Elemental  abundances  for  the .?  systems  vs  solar,  assuming  SI  =  1.0.7a  error  bars  for  the  relative  measures  are  —20%.  Diamonds — Tempel  1 ,  squares — HD 
100546,  triangles — Hale-Bopp.  The  nominal  solar  value  is  denoted  by  the  dashed  line.  The  largest  apparent  differences  are  in  the  C,  O.  and  Ca  abundances.  The  C 
abundances  may  be  affected  by  the  ISO-SWS  calibration  (see  text). 


but  metals  are  depleted.  Given  the  chemical  oddities  in  the  pho¬ 
tosphere,  it  is  not  too  surprising  that  there  might  be  some  in  the 
circumstellar  disk  as  well.  HD  has  also  been  noted  as  having  an 
unusual  mid  IR  SED  for  an  HBe  system  (Bouwman  et  al.,  2003; 
Grady  et  al.,  2005),  with  a  huge  IR  excess  due  to  abundant  small 
dust  grains.  So  the  unusual  silicate  and  carbonate  abundances 
might  be  due  to  an  unusual  starting  composition,  or  to  the  for¬ 
tuitous  circumstance  of  observing  a  system  during  the  relatively 
brief  active  disk  clearing  phase  after  planetary  core  formation, 
when  shock  and  collisional  processing  of  nebular  material  is 
active  (cf.  Wyatt  et  al.,  2007). 

On  the  other  hand,  comparing  HB  and  Tl,  the  two  Solar 
System  bodies,  a  clear  temporal  progression  can  be  seen.  The 
silicates  in  Tl  are  all  much  more  crystalline  than  in  HB,  and  the 
amount  of  PAHs  is  reduced  by  a  factor  of  25.  In  contrast,  HB’s 
composition  is  closer  to  the  mix  of  silicates  and  PAHs  in  the 
ISM  than  Tl .  The  difference  in  the  two  comets  cannot  be  due  to 
differences  in  the  location  of  formation,  as  HB  originated  in  the 
denser  and  warmer  giant  planet  region,  where  rates  of  process¬ 
ing  would,  if  anything,  be  higher  than  in  the  trans-neptunian 
region  where  Tl  formed,  leading  to  a  more  processed  body. 

5.3.  Composition  of  the  3  systems — materials  expected  in  the 
PSN  and  YSOs 

The  range  of  temperatures  of  formation  for  the  most  sta¬ 
ble  common  chemical  compounds  found  in  the  PSN,  comets, 
and  YSOs  is  as  follows.  In  general,  the  simple  metal  oxides, 
C  grains,  and  CO  are  the  most  stable  molecules  and  are  ob¬ 
served  in  the  atmospheres  of  asymptotic  giant  branch  (AGB) 
stars  at  temperatures  up  to  2500  K.  They  are  followed  in  sta¬ 
bility  by  the  silicates  which  condense  in  the  outer  atmospheres 


and  thick  winds  of  AGB  stars  at  temperatures  of  1000-1400  K. 
Next  in  stability  are  the  phyllosilicates  and  native  Fe/Fe-Ni 
metal,  forming  in  the  range  600-900  K.  The  carbonates  form 
at  600-700  K  and  the  sulfides  arc  stable  below  ~600  K.  Water 
ice  is  stable  at  temperatures  <200  K  in  vacuum,  while  the  ices 
of  methane,  carbon  monoxide,  and  carbon  dioxide  evaporate  at 
much  lower  temperatures,  below  100  K  in  vacuum. 

This  flowdown  of  materials  assumes,  of  course,  that  they 
formed  slowly  as  condensates  from  a  gas  where  all  pre-existing 
dust  had  been  vaporized.  Due  to  the  survival  of  many  types  of 
pre-solar  grains,  we  know  this  process  did  not  occur  entirely 
uniformly.  Dynamic  processes,  such  as  mixing,  shocks,  and  tur¬ 
bulence  can  create  local  non-homogeneities,  complicating  the 
situation.  Stratification  of  materials  and  reaction  processes  are 
also  expected,  as  the  Sun,  the  energy  source  for  the  formation 
of  materials,  was  located  at  the  center  of  the  system  rather  than 
spaced  uniformly  throughout  It  was  also  a  source  of  stochastic 
heating  by  high-energy  photons.  Despite  these  potential  com¬ 
plications,  the  condensation  sequence  is  still  a  useful  construct 
to  consider  as  a  measure  of  the  relative  stability  of  individual 
minerals. 

In  the  PSN,  the  dominant  species,  of  course,  was  H2.  Thus 
the  early  PSN  was  a  highly  reducing  environment.  Any  species 
not  highly  refractory  and  not  already  formed,  like  the  metal 
oxides  and  silicates  around  AGB  stars  and  novae,  would  be 
predominantly  in  the  form  of  CH4,  C2H6,  H2O,  H2S,  NH3, 
etc.  That  is,  simple  molecules  were  the  most  abundant,  and  the 
majority  of  species  were  compounds  of  hydrogen.  Some  excep¬ 
tions  to  this  rule  exist,  such  as  the  extremely  stable  N2,  CO  and 
CO2  molecules.  Once  formed  in  molecular  clouds,  they  can  last 
for  long  periods  of  time  in  reasonably  quiescent  (e.g.,  shock- 
free  and/or  low  UV- field)  environments.  Carbon  grains  formed 
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of  graphite  like  lattices  are  also  extremely  stable  once  formed. 
Hydrogenation  of  the  surface  can  occur,  though,  and  some  of 
this  surface  hydrogenation  will  lead  to  PAH  formation. 

There  was  little  or  no  free  O  in  the  form  of  O2  to  perform  the 
redox  reactions  so  familiar  to  terrestrial  and  planetary  geology. 
The  most  abundant  oxidant  available  was  instead  sulfur,  ob 
tained  from  the  labile  H2S  molecule.  The  formation  of  sulfides 
from  native  Fe  and  Fe-Ni  mixtures  was  an  oxidation  reaction. 
The  fact  that  no  sulfates  have  been  found  in  primitive  materi¬ 
als  reflects  the  lack  of  reactive  oxygen.  Sulfates  are  the  normal 
end  product  of  oxidative  attack  on  sulfide  materials,  but  the  ox¬ 
idative  attack  on  native  Fe  stopped  at  the  sulfide  level  in  the 
PSN.  Carbonates,  on  the  other  hand,  can  be  formed  directly 
from  CO2,  water,  and  silicates.  The  oxidation  of  the  carbon  and 
silicon  by  oxygen,  and  the  binding  of  the  oxygen  in  the  compo¬ 
nent  molecules,  has  already  occurred. 

In  sum,  the  solid  state  species  we  expect  to  see  in  bulk  in 
comets  and  YSOs  is  as  follows: 

•  C  in  carbon  grains,  organic  ices,  and  PAHs; 

•  O  in  silicates,  oxides,  and  water; 

•  S  in  sulfides  and  H2S; 

•  Si  in  silicates,  generally  in  the  form  of  olivines  and  pyrox¬ 
enes; 

•  Fe  in  silicates,  Fe-Ni  grains  and  Fe~NiS*; 

•  Mg  in  silicates,  MgS,  and  MgO; 

•  Al,  Ca  in  silicates  and  as  simple  oxides. 

This  closely  corresponds  to  the  range  of  materials  making  up 
the  compositional  model  for  the  Spitzer  observations  of  the  DI 
experiment  (Lisse  et  al.,  2006),  found  by  systematic  fitting  of 
laboratory  mineral  spectra  to  the  T1  ejecta  spectra,  and  of  the 
species  found  in  HB  and  HD  100546  in  this  work. 

5.4.  Compositional  trends  in  the  Tl/HB/HD 100546  spectra 

The  results  from  the  new  modeling  raise  a  number  of  issues. 
Because  the  sample  size  (2  comets  and  l  YSO)  is  so  small,  we 
present  each  as  a  question.  Definitive  answers  are  impossible 
for  a  small  dataset  with  known  selection  biases.  On  the  other 
hand,  examination  of  the  dataset  can  be  used  to  suggest  trends 
and  processes  that  may  be  verified  in  future  studies  with  a  larger 
numbers  of  objects. 

5.4.1.  Do  silicates  change  with  system  age? 

The  crystallinity  and  Fe/Mg  composition  balance  of  the  sil¬ 
icates  provides  a  high  temperature,  chemically  robust  record 
of  the  formation  history  of  the  cometary  material.  For  exam¬ 
ple,  pyroxenes  are  easier  to  modify  thermally  than  olivines 
(Hallenbeck  et  al.,  1998;  Nuth  and  Johnson,  2006).  From  the 
DI  experiment  the  T1  silicate  composition  is  more  than  80% 
crystalline,  has  an  olivinerpyroxene  ratio  of  about  1 ,  and  has  rel¬ 
atively  good  equilibration  between  the  Fe-  and  Mg-rich  in  both 
the  olivine  and  pyroxene  phases.  The  silicates  in  HD  are  a  very 
simple  system — only  Mg-rich  crystalline  olivine,  amorphous 
pyroxene,  and  some  crystalline  pyroxenes  are  present,  with 
about  50%  total  crystalline  material  and  an  olivinerpyroxene  ra¬ 


tio  of  about  1.  The  HB  system  seems  intermediate  between  the 
two,  with  ~67%  crystalline  olivines  and  pyroxenes,  relatively 
more  Fe/Ca/Mg-equilibiated  pyroxenes,  a  goodly  amount  of 
amorphous  olivine,  but  again  only  Mg-rich  crystalline  olivine. 
The  silicate  evidence  implies  that  the  dust  in  T1  has  been  most 
heavily  processed  of  the  three  systems. 

To  test  if  the  different  apparent  nature  of  the  T1  silicates 
is  due  to  the  different  spectral  resolutions  and  sensitivities  for 
the  SWS  and  IRS  instruments,  we  smoothed  the  T1  spectrum  to 
match  the  “rounder*  features  of  the  ISO  HB  data.  We  recovered 
the  same  mix  of  materials  as  found  for  the  unsmoothed  data, 
including  the  Fe-rich  silicates  not  found  in  HB.  Also,  the  pre¬ 
liminary  results  of  Stardust  for  Comet  Wild  2,  a  Jupiter  family 
comet  (JFQ  like  Tl,  confirm  the  presence  of  Mg-rich  olivines 
and  Ca-  and  Fe-rich  pyroxenes  in  the  Wild  2  dust  (Brownlee  et 
al.,  2006b). 

If  silicates  were  processed  by  annealing  near  the  Sun,  then 
this  implies  an  efficient  mechanism  for  transferring  silicaceous 
dust  from  the  inner  Solar  System  to  the  K-E  belt  for  some 
time  after  the  Oort  cloud  comets  such  as  HB  formed  in  the 
giant  planet  region  (Gail,  2004;  Keller  and  Gail,  2004).  How¬ 
ever,  efficient  radial  transport  from  the  inner  nebula  should  end 
at  about  the  time  of  giant  planet  core  formation,  at  ~10  Myr 
(Boss,  2004).  The  density  of  material  in  the  giant  planet  region 
was  much  higher  than  at  the  edges  of  the  PSN  in  the  K-E  belt 
where  JFCs  formed.  Thus,  Oort  cloud  comets  like  HB,  on  the 
whole,  probably  formed  before  JFCs  like  Tl.  We  can  “date” 
the  average  Oort  cloud  comet  at  less  than  10  Myr,  the  time  es¬ 
timated  for  the  giant  planet  core  formation — the  giant  planet 
planetesimals  were  formed  first,  before  the  giant  planets  them 
selves.  This  leaves  the  relatively  short  time  between  formation 
of  the  planetesimals  and  the  first  planet  cores  as  the  source  of 
the  difference  in  the  silicates  in  Tl  and  HB,  even  when  we  also 
consider  nebular  shock  processing  and  evolutionary  processing 
due  to  repeated  solar  passages. 

The  shock  theory  of  Harker  and  Desch  (2002)  argues  that 
crystalline  silicates  should  start  forming  at  the  time  of  giant 
planet  core  creation  and  growth.  Shocks  had  only  until  the  time 
of  orbital  ejection  of  the  planetesimal  into  the  Oort  cloud  to 
work  on  the  material  of  HB.  The  material  in  the  slower  forming 
Tl  could  have  been  processed  until  the  giant  planets  finished 
accreting.  The  existence  of  the  primitive  chondrules,  which  ap¬ 
pear  to  have  been  formed  from  liquid  phase  solidification  of 
molten  PSN  material,  provides  strong  evidence  for  shocks  in 
the  PSN.  If  shocks  formed  the  primitive  chondrules  and  altered 
the  silicaceous  material  in  Tl,  then  the  mineral  composition  of 
JFCs  like  Tl  should  be  closer  to  chondtitic  than  Oort  cloud 
comets  like  HB  which  formed  in  the  pre-shock  era.  However, 
while  there  is  good  agreement  between  the  Tl  and  Cl  chon- 
dritic  abundances,  there  is  equally  good  agreement  between  the 
HB  and  the  Cl  abundances.  The  low  Fe/Mg  crystalline  olivine 
ratio,  low  particle  strength,  and  high  porosity  also  argue  for 
mainly  thermal  processing  and  against  much  shock  process¬ 
ing  of  the  material  in  Tl,  as  suggested  by  Lisse  et  al.  (2006). 
On  the  other  hand,  as  noted  above,  shocks  may  very  well  be 
strongly  processing  the  HD  silicate  material,  reforming  it  into  a 
mix  of  Mg-rich  olivine  and  amorphous  pyroxene.  To  have  left 
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the  other  materials  in  the  system  with  a  comet-like  composi¬ 
tion  while  altering  the  highly  refractory  silicates  takes  some 
special  conditions.  It  would  seem  to  require  destruction  of  the 
materials  in  the  inner  disk  cavity  (r*  <  13  AIJ)  and  efficient 
transport  of  the  shock  processed  material  into  the  edge  of  the 
disk  at  r*  ~  13  AU,  adding  to  and  overwhelming  the  local,  un¬ 
processed,  “native”  silicate  composition. 

The  Mg/Fe  silicate  equilibration  seen  in  T1  could  have  oc¬ 
curred  since  cometary  formation,  in  the  repeated  close  pas¬ 
sages  by  the  Sun,  building  up  a  mantle  of  material  substantially 
different  from  the  bulk  interior.  This  process  would  also  ap¬ 
pear  to  have  caused  an  equilibration  of  the  Mg/Fe  abundance 
in  the  sulfide  phases  as  well  (Section  5.4.4).  However,  it  is 
unlikely  that  Fe/Mg  cation  migration  could  occur  at  the  low 
temperatures  found  in  cometary  bodies.  It  clearly  has  not  in 
HB,  while  chemistry  leading  to  carbonate  and  phyllosilicate 
formation  has.  Further,  there  is  now  abundant  evidence  from 
the  DI  mission  that  the  average  material  on  T1  was  thermally 
processed  to  depths  of  less  than  1  m  (Groussin  et  al.,  2007; 
A' Hearn  et  al.,  2007)  while  the  DI  impactor  ejected  material 
from  depths  of  10-30  m  deep.  Barring  any  extensive  reworking 
of  the  interior  material  by  impact  processing,  we  thus  conclude 
that  there  was  an  exceptional  amount  of  processing  of  nebular 
material  in  the  PSN  between  the  time  of  HB’s  formation  and 
the  end  of  the  era  of  giant  planet  formation,  when  substantive 
inputs  to  the  T1  composition  would  have  ended. 

5.4.2.  Are  PAHs  as  common  in  nebulae  and  comets  as  in  the 
ISM? 

PAHs,  while  not  the  major  carriers  of  carbon  in  the  Solar 
System  and  ISM  (they  are  observable  mainly  due  to  their  pro¬ 
nounced  UV-driven  fluorescence)  are  nevertheless  important 
species  chemically  and  useful  tracers  of  processing  in  stellar 
environments.  In  the  diffuse  ISM,  the  abundance  of  PAHs  vs  sil¬ 
icates  is  1-10%,  orders  of  magnitudes  higher  than  what  is  found 
in  comets  (Li  and  Draine,  2001 ;  Li  and  Draine,  2007.  On  PAHs 
is  comets,  in  preparation).  This  suggests  that  while  PAHs  could 
be  part  of  the  earliest  material  input  into  the  PSN,  they  might 
have  been  removed  or  destroyed  as  the  system  evolves.  Does 
the  removal  of  PAHs  begin  as  a  molecular  cloud  core  collapses 
into  the  protostellar  disk,  or  does  it  happen  later,  as  the  central 
source  turns  on  and  irradiates  the  circumstellar  material?  The 
strong  PAH  signal  in  HD  and  many  other  YSOs  (Sloan  et  al., 
2005)  argues  that  PAHs  often  survive  infall  and  are  destroyed 
later  (e.g.,  photodissociated  in  the  PSN  after  they  emerge  out  of 
the  sublimating  ice  mantles  of  dust  grains). 

The  fact  that  other  parts  of  the  HD  spectrum  seem  to 
match  cometary  spectra  well  has  been  used  in  the  past  to  ar¬ 
gue  that  PAHs  should  be  common  in  comets.  As  mentioned 
in  the  HB  compositional  modeling  section,  the  comet  was 
found  to  be  rich  in  simple  organic  species  (Biver  et  al.,  1999; 
Bockelee-Morvan  et  al.,  2000),  so  it  is  not  surprising  that  car¬ 
bonaceous  PAHs  were  found  using  the  detailed  modeling  de¬ 
scribed  in  this  work.  Ionized  PAHs  are  found  in  all  3  systems 
in  this  study,  and  are  beginning  to  seem  ubiquitous  at  the  thou¬ 
sand  ppm  level,  as  the  mid  IR  spectra  from  two  other  comets 
observed  by  Spitzer  also  show  evidence  for  PAHs  (Lisse  et  al.. 


2007,  in  preparation).  The  markedly  stronger  emission  features 
would  seem  to  aigue  for  much  more  PAH  abundance  in  HD 
than  HB,  but  the  abundances  in  Tables  2  and  3  indicate  they 
are  at  roughly  the  same  abundance  level,  2  x  10  2  molar  ratio. 
This  seeming  contradiction  can  be  resolved  when  the  mecha¬ 
nism  driving  the  PAH  emission,  UV  pumping  by  the  G2  V  Sun 
or  the  Be9  V  HD  primary,  rather  than  direct  heating  by  stellar 
insolation,  is  accounted  for.  It  is  also  important  to  note  that  the 
abundance  level  of  PAHs  in  HD  and  HB  is  about  the  same  as 
in  the  ISM,  implying  that  PAHs  can  indeed  survive  infall  and 
incorporation  into  young  Solar  System  bodies. 

The  larger  abundance  of  PAHs  in  HB  compared  to  Tl,  about 
a  factor  of  25  (Lisse  et  al.,  2006),  appears  to  be  a  real  effect, 
which  we  ascribe  to  the  later  time  of  formation  of  Tl,  as  PAHs 
become  gradually  more  scarce  in  the  PSN  as  they  are  destroyed 
over  time.  As  for  the  case  of  silicate  processing,  we  believe  we 
can  rule  out  the  effects  of  different  formation  locations.  In  the 
mixing  scenario,  PAHs  would  be  destroyed  more  quickly  in  the 
giant  planet  region,  where  solar  UV  and  molecular  H2  were 
much  more  abundant,  than  in  the  K-E  belt  The  PAHs  in  the 
giant  planet  region  could  easily  have  continued  on  the  path  of 
reduction  to  still  smaller,  saturated  organics  by  photochemistry 
in  the  PSN.  This  process  would  have  taken  longer  to  perform  in 
the  relatively  dark  and  diffuse  K-E  belt  allowing  substantially 
more  PAHs  to  survive  intact  there.  The  shock  model  would  also 
argue  for  a  more  rapid  destruction  of  PAHs  in  the  dense  PSN 
material  of  the  giant  planet  region,  which  is  capable  of  support¬ 
ing  stronger  shocks  than  the  sparse  material  in  the  K-E  belt. 
The  fact  that  we  see  a  higher  relative  abundance  of  PAHs  in  HB 
than  in  Tl  argues  for  an  earlier  age  of  formation  for  HB,  as  does 
the  more  primitive  nature  of  the  silicates. 

5.4.3.  Are  carbonates  and  phyllosilicates  present? 

The  presence  of  the  minor  carbonate  and  phyllosilicate 
species  in  YSOs  or  comets  is  an  issue  currently  in  much  de¬ 
bate.  Ceccarelli  et  al.  (2002)  and  Chiavassa  et  al.  (2005)  found 
that  half  of  the  YSOs  in  their  ISO  Long  Wavelength  Spectrom¬ 
eter  (LWS;  Gry  et  al.,  2003)  study  have  carbonate  emission. 
Recent  work  by  Toppani  et  al.  (2005)  demonstrated  that  carbon¬ 
ates  can  form  from  gas-phase  silica,  water,  and  carbon  dioxide, 
although  at  relatively  high  partial  pressures  versus  those  found 
in  the  PSN.  Re-analysis  of  the  IDP  record  has  shown  evidence 
for  carbonate  materials  (Flynn  et  al.,  2005).  Comets  are  known 
to  contain  abundant  water  ice,  carbon  dioxide,  and  silicates,  all 
which  can  react  in  the  parent  body  to  form  carbonates  and  phyl¬ 
losilicates  in  planetary  environments.  Also,  there  are  no  other 
likely,  highly  abundant  sources  of  emission  in  the  wavelength 
range  6.5-7.5  pm,  where  we  find  the  carbonate  peak. 

On  the  other  hand,  the  canonical  wisdom  holds  that  carbon¬ 
ates  and  phyllosilicates  only  form  in  the  presence  of  liquid  wa¬ 
ter.  This  would  mean  that  liquid  water  must  have  been  present  in 
Tl  and  HB  at  some  point,  a  finding  surprising  to  many,  although 
new  work  by  Merk  and  Prialnik  (2006)  argues  for  abundant  liq¬ 
uid  water  in  KBOs  during  the  era  of  formation  due  to  heating 
by  Al26.  Also,  there  is  some  controversy  over  the  claim  of  car¬ 
bonate  detection  by  Ceccarelli  et  al.  (2002)  and  Chiavassa  et  al. 
(2005).  The  emission  features  used  to  determine  the  presence 
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and  amount  of  carbonates  found  in  LWS  spectra  at  90- 1 10  um 
are  relatively  broad  and  shallow.  They  depend  highly  on  the 
continuum  subtraction  in  many  cases,  and  the  results  could  be 
confused  with  emission  lines  due  to  silicates  and  phyllosilicates 
(e.g.,  Bouwman  et  al.,  2003). 

Definitive  evidence  for  carbonates  or  phyllosilicates  has  yet 
to  be  found  in  the  Stardust  samples  studied  to  date,  which  is 
worrisome.  It  is  not  clear  at  the  time  of  this  writing  how  much 
small  number  statistics,  hot  sample  collection  into  aerogel,  col¬ 
lection  of  mainly  surface  mantle  (and  thereby  altered)  material, 
or  samples  from  a  substantially  younger  comet  affect  the  Star¬ 
dust  samples  vs  the  Deep  Impact  experiment  results;  e.g.,  so  far 
only  the  most  refractory  of  species  (silicates  and  sulfides)  have 
been  recovered  intact  in  the  terminal  particles  of  the  Stardust 
tracks;  UV  photodecomposition  could  transform  carbonates  in 
the  surface  layers  of  Comet  Wild  2  into  metal  oxides  and  CO2 
(as  occurs  on  Mars,  where  no  carbonates  are  found,  despite  the 
prevalence  of  species  and  structures  formed  by  aqueous  water 
alteration  of  the  surface  material,  similar  to  where  carbonates 
are  found  terrestrially);  on  as  no  obvious  impact  craters  were 
found  on  the  surface  of  Wild  2  (Brownlee  et  al.,  2004),  the  hulk 
material  of  Comet  Wild  2  may  have  undergone  no  substantial 
impact  processing.  With  respect  to  the  latter  point,  it  should  be 
kept  in  mind  that  Deep  Impact  sampled  material  from  one  place 
at  depth  on  T1 ,  located  between  two  300  m  wide  impact  craters, 
where  reworking  of  the  comet  has  clearly  occurred.  It  is  quite 
possible  that  the  T1  ejecta,  while  definitively  from  the  interior 
of  the  comet,  contains  some  additional  species,  along  with  the 
materials  found  in  a  young  comet,  due  to  impact  driven  chemi 
cal  processing  (Lisse  et  al.,  2006). 

Given  this  uncertain  context,  we  find  it  telling  that  our  model 
detects  carbonates  in  the  highly  active  and  water  rich  HB  sys¬ 
tem  and  in  T1  as  well  as  in  the  HD  YSO  system,  but  only  after 
we  carefully  removed  the  photospheric  contribution.  This  gives 
us  confidence  that  our  model  does  not  detect  ‘false  positive*  car¬ 
bonates.  One  other  point  can  be  made  from  our  modeling:  the 
relationship  (or  lack  thereof)  between  the  parent  silicates  and 
daughter  carbonate  products.  In  HB,  the  carbonates  are  slightly 
Fe-rich,  with  a  6:5  ratio  of  siderite  to  magnesite,  with  about 
15%  abundance  by  surface  area.  For  Tl,  the  ratio  is  5:3,  with 
about  5%  abundance.  Both  ratios  mimic  the  pyroxene,  but  not 
the  olivine,  Fe:Mg  ratios.  For  HD  we  find  only  Mg-rich  car¬ 
bonate  at  about  3%  abundance,  which  agrees  with  the  large 
preponderance  of  Mg-rich  olivine,  but  not  the  relatively  Fe-rich 
pyroxenes.  The  unusually  large  amount  of  amorphous  pyrox¬ 
ene  in  HD  may  be  a  clue  to  unusual  pyroxene  chemistry  in  this 
system.  Also,  the  effect  of  shocks  and  collisions  on  carbonate 
materials  is  unclear,  although  we  would  expect  them  to  be  rela¬ 
tively  “unrefractory”  vs  the  other  species  in  our  compositional 
mix,  decomposing  to  a  mix  of  carbon  dioxide  and  metal  oxide. 
More  detailed  IR  observations  of  comets  and  YSO  systems  are 
clearly  needed  to  answers  these  questions. 

5.4.4.  What  can  he  learned  from  the  sulfide  mineralogy  in  the 
3  systems? 

Forming  at  ~600  K,  sulfides  have  the  promise  of  being  ex¬ 
tremely  useful  diagnostics  of  the  local  PSN  environment  (Keller 


et  al.,  2002).  The  major  difference  among  the  three  systems  is 
that  HB  and  HD  100546  are  much  better  fit  by  an  Fe-rich  sul¬ 
fide  component  with  stoichiometry  Feo.9MgO.1S  as  opposed  to 
the  Feo.5Mgo.5S  found  in  the  Tl  ejecta.  From  this  and  the  deter¬ 
mination  that  HB  and  HD  100546  are  less  chemically  evolved 
than  Tl,  we  infer  that  the  less  equilibrated  sulfides  are  more 
primitive.  Preliminary  evidence  from  the  Stardust  mission  is 
that  the  Comet  81P/Wild  2  sulfides  occur  in  the  form  of  non- 
stnichiometric  FeSx  (with  admixtures  of  Ni;  Keller  et  al.,  2006; 
Zolensky  et  al.,  2006).  This  mineral  system  is  currendy  poorly 
understood  by  ground  based  laboratory  studies  and  as  a  con¬ 
sequence  is  not  well  supported  in  our  model,  but  does  appear 
to  be  more  closely  related  to  the  Fe-rich  ningerites  found  for 
HB  and  HD  in  this  work.  In  addition,  some  quantity  of  S  bear 
ing  molecules  such  as  H2S  or  SO2  may  not  react  with  available 
metals  and  thus  may  not  become  incorporated  in  the  dusty/icy 
aggregates  that  accrete  to  form  the  comets.  Until  better  labora¬ 
tory  and  Stardust  measurements  are  made,  the  sulfide  question 
will  remain  open. 

6.  Conclusions 

1.  The  same  model  used  to  fit  the  Tl  ejecta  composition 
fit  the  HB  and  HD  systems  well.  The  model  is  based  as 
much  as  possible  on  experimental  observation  and  materi¬ 
als  analysis  of  sample  returns.  Dominant  features  for  the 
systems  found  in  other  models,  like  the  dominance  of  Mg- 
rich  olivines  and  water  ice  in  HB  and  the  presence  of  PAHs 
and  Mg-rich  olivines  in  HD,  were  reproduced. 

2.  Particle  temperature  and  grain  size  effects  must  be  included 
in  the  models  to  accurately  determine  compositional  min¬ 
eralogy  with  mid-IR  spectra.  The  material  in  Tl  yielded 
very  different  emission  spectra  when  it  was  emitted  natu 
rally  in  large  aggregates  (ambient  coma)  or  separated  into 
a  fine  powder  by  the  DI  experiment  (ejecta). 

3.  The  detailed  mid-IR  spectra  of  the  Tl  ejecta,  the  Comet  HB 
coma,  and  the  HD  nebula  are  qualitatively  similar.  How¬ 
ever,  there  are  important  differences  between  the  systems 
as  well  which  allow  for  a  comparative  mineralogy  study. 

4.  The  material  incorporated  in  HB  appears  to  be  less  proc¬ 
essed  than  that  in  Tl ,  as  determined  from  the  differences  in 
their  crystalline  silicate  content  The  relative  abundance  of 
PAHs  in  HB  and  in  Tl  can  be  explained  by  a  net  destruc 
tion  over  time  of  PAHs  in  the  PSN,  along  with  an  earlier 
formation  time  for  comets  in  the  giant  planet  region  vs  the 
K-E  belt. 

5.  Both  comets  contain  carbonates  and  phyllosilicates,  as  does 
HD 100546.  The  cometary  carbonate  Fe:Mg  composition 
follows  the  pyroxene  composition  in  the  comet  systems. 

6.  The  HD  SED  can  be  explained  by  emission  dominated  by 
a  single  mix  of  material  co- located  at  ~13  AU.  The  two 
different  apparent  temperature  regimes  are  due  to  the  pres 
ence  of  abundant  hot  refractory  dust  mixed  with  abundant 
sublimation  cooled  water  ice  particles. 

7.  The  material  composition  of  HD  100546  seems  to  be  very 
similar  to  that  of  the  two  comet  systems,  with  the  exception 
of  the  absence  of  amorphous  olivine  and  the  superabun- 
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dance  of  amorphous  pyroxene.  The  presence  of  abundant 
PAHs  and  water  ice  implies  that  there  is  little  processing 
of  the  nebular  material.  The  unusual  silicate  composition 
may  be  due  to  either  a  different  starting  composition  for 
the  HD  nebula  vs  the  PSN,  or  due  to  shock  and  collisional 
processing  of  material  inside  13  AU  caused  by  due  to  planet 
formation  and  disk  clearing,  followed  by  transport  into  the 
inner  edge  of  the  surrounding  disk. 

8.  PAHs  are  found  in  all  3  systems  studied.  The  abundance  of 
PAHs  in  HB  and  HD  is  roughly  the  same,  ~2%  vs  the  bulk, 
similar  to  the  level  seen  in  the  ISM.  T1  is  depleted  in  PAHs, 
by  a  factor  of  25. 
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